The source of the "third flow component" 
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A model for energy, pressure and flow velocity distributions at the beginning of 
relativistic heavy ion collisions has been developed to be used as initial condition 
for hydrodynamical calculations. The results show that at the early stages QGP 
forms a tilted disk, such that the direction of the largest pressure gradient stays 
in the reaction plane, but deviates from both the beam and the usual transverse 
flow directions. Such initial conditions may lead to the creation of the "antiflow" 
or the "third flow component" [jj. 



1 Introduction 

Fluid dynamics is probably the most frequently used model to describe heavy 
ion collisions. It assumes local equilibrium, i.e. the existence of an Equation of 
State (EoS) , relatively short range interactions and conservation of energy and 
momentum as well as of conserved charge(s). Thus, in ultra relativistic colli- 
sions of large heavy ions, especially if Quark-Gluon Plasma (QGP) is formed, 
one-fluid dynamics is a valid and good description for the intermediate stages of 
the reaction. Here interactions are strong and frequent, so that other models, 
(e.g. transport models, string models, etc., assuming binary collisions, with 
free propagation of constituents between collisions) have limited validity. On 
the other hand, the initial and final, Freeze-Out (FO), stages of the reaction 
are outside the domain of applicability of the fluid dynamical model. 

We believe that the realistic and detailed description of an energetic heavy 
ion reaction requires a Multi Module Model, where the different stages of the 
reaction are each described with suitable theoretical approach. It is important 
that these Modules are coupled to each other correctly: on the interface, which 
is a 3 dimensional hyper-surface in space-time, all conservation laws should be 
satisfied and entropy should not decrease. These matching conditions were 
worked out and studied for the matching at FO in detail in refs. 
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The final FO stages of the reaction, after hadronization, can be described 
well with kinetic models, since the matter is already dilute. 

The initial stages are more problematic. Frequently two or three fluid 
models are used to remedy the difficulties and to model the process of QGP 
formation and thermalization&lm Here the problem is transferred to the de- 
termination of drag-, friction- and transfer- terms among the fluid components, 
and a new problem is introduced with the (unjustified) use of EoS in each com- 
ponent in noncquilibrated situations, where EoS does not exist. Strictly speak- 
ing this approach can only be justified for mixtures of noninteracting ideal gas 
components. Similarly, the use of transport theoretical approaches assuming 
dilute gases with binary interactions is questionable, since, due to the extreme 
Lorentz contraction in the CM. frame, enormous particle and energy densi- 
ties with the immediate formation of perturbative vacuum should be handled. 
Even in most parton cascade models these initial stages of the dynamics are 
just assumed in form of some initial condition with little justification behind. 

Our goal in the present work is to construct a model, based on the re- 
cent experience gained in string Monte Carlo models and in parton cascades, 
describing energy, pressure and flow velocity distributions at the beginning of 
rclativistic heavy ion collisions. 



2 The effective string rope model 

One of the important results of the last decadeis that all string-models had to 

a' ltroduce new, energetic objects: string ropesEra, quark clustersQ, fused strings 
, in order to describe the abundant formation of massive particles like strange 
antibaryons. Based on this, we describe the initial stages of the reaction hi 
the framework of classical (or coherent) Yang-Mills theory, generalizing ref. E3 
assuming larger field strength (string tension) than in ordinary hadron-hadron 

Slhsions. In addition we now satisfy all conservation laws exactly, while in ref. 
infinite projectile energy was assumed, and so, overall energy and momentum 
conservation was irrelevant. 

First of all, we would create a grid in [x, y] plane (z - is the beam axes, [z, x] 
- is the reaction plane) and describe the nucleus- nucleus collision in terms of 
steak-by-streak collisions, corresponding to the same transverse coordinates, 
{xi,yj}. We assume that baryon recoil for both target and projectile arise 
from the acceleration of partons in an effective field F^" , produced in the in- 
teraction. Of course, the physical picture behind this model should be based on 
chromoelectric flux tube or string models, but for our purpose we consider F^" 
as an effective abelian field. A single phenomenological parameter describing 
this field must be fixed from comparison with experimental data. For detailed 
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description we would call the attention of the reader to refs. 




Briefly speaking, after two streaks cross each other they create a "string 
rope" with constant string tension a in the middle. The typical values of the 
string tension, a, are of the order of 10 GeV/fm, and these may be treated 
as several parallel strings or " string rope" . In further evolution this " string 
rope" accommodates more and more energy by stopping the initial target and 
projectile streaks. If we let system evolve, the "string rope" will stop streaks 
completely, then will accelerate them backward, and the analog of Yo-Yo mo- 
tion may occur. We do not solve simultaneously the kinetic problem leading 
to parton equilibration, just assume that the arising friction is such that the 
heavy ion system will be an overdamped oscillator, i.e. yo-yoing of the two 
heavy ions will not occur. This assumption is based on recent string and parton 
cascade results. 

To proceed one has to make more assumptions. We assume that the final 
result of collisions of two streaks, after stopping the string's expansion and after 
its decay, is one streak of the length Alf with homogeneous energy density 
distribution, e/, and baryon charge distribution, pf, moving like one object 
with rapidity yj. We assume that this is due to string-string interactions and 
string decays. The homogeneous distributions are the simplest assumptions, 
which may be modified based on experimental data. Its advantage is a simple 
expression for e/ , pf, y/. 

The final energy density, baryon density and rapidity, e/, pf and j//, 
should be determined from conservation laws. Unfortunately, the assump- 
tions we made above oversimplify the real situations and do not allow us to 
satisfy exactly all the conservation lawsOilj. The reason for this is well known 
and has been discussed in the Refs. Era: two possible definitions of the flow, 
Eckart's and Landau's definition. If we are following the energy flow, we sat- 
isfy exactly the energy and momentum conservation and violate the baryon 
current conservation. Otherwise if we are drifted by baryon flow, we violate 
the energy-momejatrrjiLconservation. 

In the refs. Eirli3ll3 we choose the Landau's definition of the yf, i.e. the 
exact conservation of the energy and momentum. 



3 Initial state of QGP 

Let us review the results of our calculations. We are interested in the shape 
of QGP formed, when the final streak formation is finished (at least in the 
most central regionliil) and their matter is locally equilibrated. This will be 
the initial state for further hydrodynamical calculations. The time, where we 
expect a local equilibrium, is also a parameter of our model - it should be large 
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Au+Au, e n = 100 GeV/nucI, b = 0.5^ +rA A=0.051 
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Figure 1: The Au+Au colKaifMujqeo = 100 GeV/nucl, b = 0.5(Ki + R 2 ), A = 0.051 (param- 
eter A introduced in refs. PHri ' M to determine the "string rope" tension), y = (ZX plane 
through the centers of nuclei). We would like to notice that final shape of QGP volume is 
a tilted disk si 45°, and the direction of the fastest expansion will deviate from both the 
beam axis and the usual rfcsansverse flow direction, and might be a reason for the third flow 
component, as argued inPl. 



enough to allow the final streak formation, but we can't wait too long, since 
the transverse expansion can't be neglected in this case. 

We see in Figs. 0, | that QGP forms a tilted disk for 6/0. So, the di- 
rection of fastest expansion, the same as largest pressure gradient, will deviate 
from both the beam axis and the usual transverse flow direction and the new 
flow component, called "antiflow" or "third flow component", may appear in 
the reaction plane. With increasing beam energy the usual transverse flow is 
getting weaker, while this new flow component is strengthened. The mutual 
effect of the usual directed transverse flow and this new "antiflow" or "third 
flow component" leads to an enhanced emission in the reaction plane. This 
was actually observed and widely studied earlier and was referred to as " elliptic 
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flow" . 



Au+Au, e n = 100 GeV/nucI, b = 0.25^ +r A A=0.051 




Figure 2: The same as Fig. ^, but b = 0.25(ri + T2)- We see that for more central collisions 
the energy density is much largpn The QGP volume has a shape of tilted disk and may 
produce a third flow component EJ. 



4 Conclusions 

Based on earlier Coherent Yang-Mills field theoretical models and introducing 
effective parameters based on Monte-Carlo string cascade and parton cascade 
model results, a simple model is introduced to describe the pre fluid dynamical 
stages of heavy ion collisions at the highest SPS energies and above. 

Contrary to earlier expectations, — based on standard string tensions of 
1 GeV/ fm which lead to the Bjorken model type of initial state, — the effective 
string tension introduced for collisions of massive heavy ions, as a consequence 
of collective effects related to QGP formation, appears to be of the order of 
10 GeV/ fm and consequently causes much less transparency. The resulting 
initial locally equilibrated state of matter in semi central collisions takes a 
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rather unusual form, which can be then identified by the asymmetry of the 
caused collective flow. 

Detailed fluid dynamical calculations as well as flow experiments at semi 
central impact parameters for massive heavy ions arc needed at SPS and RH1C 
energies to connect the predicted special initial state with observables. 
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